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I. INTRODUCTION 



Lattice QCD has made impressive progress. As fully dynamical simulations for a wide 
range of observables are carried out at lighter quark masses, our confidence in lattice QCD 
as a predictive tool grows. There are a number of fermion discretizations used on the lattice. 
These are at various stages of development, and are fraught with disparate difficulties. Fur- 
ther confidence in lattice QCD will eventually be built through comparisons of simulations 
employing different lattice fermions. 

Presently lattice QCD studies using dynamical staggered fermions [l[ reach smaller quark 
masses as compared with other lattice fermions. The publicly available MILC configura- 
tions 0, moreover, have launched dynamical staggered fermions as readily accessible for 
lattice calculations [E HI IE IE 0, HI- The fourth-root trick is employed by these simula- 
tions to reduce the number of so-called taste degrees of freedom. In the continuum limit, 
there are no taste-changing interactions and the trick is kosher. On the lattice at finite a, 
however, the question of locality emerges and the trick remains controversial, for various 
recent investigations see 0, HO, U, IE IE IE IE IE 12 1- Nonetheless, even with improved 
staggered actions the discretization effects are surprisingly large. Chiral perturbation theory 
has been extended to staggered actions to control the systematic errors associated with the 
continuum extrapolation IE IE, IE 21, IE HI- 

Simulations that employ Ginsparg- Wilson fermions 24]|, on the other hand, do not suffer 
from the precarious theoretical situation of staggered fermions. 1 These fermions, more- 
over, have an exact chiral symmetry (which in practice is limited by how well the overlap 
operator 2E IE 22j is approximated, or how well the domain wall fermion jiE IE IE) is re- 
alized). Computationally, however, Ginsparg- Wilson fermions are numerically quite costly. 
Mid-ground between these two lattice fermions can be found. Ginsparg- Wilson quarks can 
be utilized for the valence quarks, where they are computationally less demanding. One can 
then calculate correlation functions in the background of the existing staggered sea of the 
MILC configurations. This efficacious solution has costs comparable to quenched Ginsparg- 
Wilson simulations, and numerical investigations in mixed-action lattice QCD have been 
undertaken [3]], I^E I^E 0|- Mixed-action chiral perturbation theory has recently been 
formulated in the meson sector [35j] to analyze the lattice-spacing dependence of meson 
observables. 

In this work, we construct the baryon chiral Lagrangian for a mixed lattice action consist- 
ing of Ginsparg- Wilson valence quarks and staggered sea quarks. We apply this Lagrangian 
to the computation of octet baryon masses and magnetic moments. To address the effects 
of finite lattice spacing, one formulates the underlying lattice theory and matches it onto a 
chiral effective theory j36|. To do so, we utilize a dual expansion in the quark masses and 
lattice spacing. We assume a hierarchy of energy scales 



m q < Aqcd < -, 
a 



(1) 



In this work, we assume the validity of the fourth-root trick. Comparing lattice data with the behavior 
of observables predicted from staggered chiral perturbation theory could yield empirical evidence for or 
against the trick. 
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and further choose the power counting scheme 



£ 2 ^ I m q /A QC D ^ 

a ^QCD 



which is currently relevant for simulations employing improved staggered quarks [6(. The 
resulting quark mass and lattice spacing dependent expressions will be useful for the analysis 
of numerical results from simulations of baryon observables in mixed-action lattice QCD. 

Our presentation has the following organization. First we briefly review the Symanzik 
Lagrangian for the mixed lattice action in Section [H] Here we recall the form of the chiral 
Lagrangian in the meson sector [35|, and cite the relevant details for our calculation. The 
baryons are then included in the theory. In Section ITTT1 we determine the magnetic moments 
to 0(e) in the combined expansion; while in Section HVl we calculate the masses of the octet 
baryons up to 0(e 3 ) . These calculations at their respective orders include the leading non- 
analytic dependence on the quark masses, and are shown to be devoid of taste-symmetry 
breaking. Taste-symmetry breaking interactions contribute at 0(e A ) from loops that scale 
generically as 0(a 2 m q \ogm q ) and 0{m 2 \ogm q ). A summary ends the paper ( Section IVjl. 
For completeness we include a discussion of the finite volume corrections in the Appendix. 



II. CHIRAL LAGRANGIAN 

Before including baryons into the mixed-action chiral Lagrangian for Ginsparg- Wilson 
valence quarks and staggered sea quarks, we first give a brief review of the Symanzik La- 
grangian. Next we recall the form of the chiral effective theory in the meson sector, and 
list results pertinent for the calculations in this paper. The relevant pieces of the baryon 
Lagrangian are then detailed. 



A. Symanzik Lagrangian 

The lattice action can be described in terms of a continuum effective field theory. This 
effective theory is described by the Symanzik action (37, 38[, which is built from continuum 



operators and is based on the symmetries of the underlying lattice theory. The Symanzik 
Lagrangian is organized in powers of the lattice spacing a, namely 

£ Sy m = C + a£ (5) + a 2 £ (6) + a 3 £ (7) + . . . , (3) 

where C^ 1 ' represents the contribution from dimension-n operators. 2 The symmetries of the 
mixed lattice action are respected by the Symanzik Lagrangian £sym order-by-order in a. 
In the continuum limit, a — > 0, only the operators of L survive. We consider here the case 
of a mixed action in partially quenched QCD (PQQCD). This type of action allows not 
only for the valence and sea quarks to have different masses, but to be different types of 
lattice fermions. An important feature of mixed action theories concerns the general lack 
of symmetry between the valence and sea sectors [3^. Because different types of lattice 



2 One should note that not all a-dependence is parametrized in Eq. J2J. The coefficients of operators in 
D- n > depend on the gauge coupling, and hence can have a weak logarithmic dependence on a. 
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fermions are used in each sector of the theory, the flavor symmetry of the mixed action is 
generally a direct product of the separate symmetries in the valence and sea sectors. 

The underlying lattice action we consider is built from three flavors of Ginsparg- Wilson 
valence quarks and three flavors of staggered sea quarks. In the continuum limit, the La- 
grangian £ is just the partially quenched action, 3 namely 

C = QIpQ + Qm q Q, (4) 

where the quark fields appear in the vector 

Q = (Uid,sJ 1 J 2 j3,j4ih,l2,h,h,r 1 ,r 2 ,r 3 ,r 4 ,u,d,s) T . (5) 

Notice that fermion doubling has produced four tastes for each flavor (j,l,r) of staggered 
quark. In a partially quenched generalization of the isospin limit, the mass matrix is given 
by 

m q = diag(m u , m u , m s , m^, m^, m r f j, m u , m u , m s ), (6) 

with £/ as the 4x4 taste identity matrix. In the massless limit, the Lagrangian C has a 
graded chiral symmetry of the form S77(15|3)l ® SU(15\3)r. For our discussion below, it is 
useful to define projection operators for the valence (V) and sea (S) sectors of the theory: 
V v = diag(l, 1, 1, 0, . . . , 0, 1, 1, 1), and V s = (0, 0, 0, lfr, If/, 0, 0, 0). 

The Ginsparg- Wilson sector of the theory possesses an exact chiral symmetry in the limit 
of zero quark mass j^. Thus there can be no operators in involving only valence quarks 
because the only dimension-5 operator (after field redefinitions [43]) is a chiral symmetry 
breaking quark bilinear. Also there are no dimension-5 operators built from just staggered 
quark fields |44j, |45j . Finally the symmetries of the mixed action forbid bilinear operators 
formed from one valence and one sea quark. Thus = 0. 

Next we consider the dimension-6 operators, of which there are purely gluonic opera- 
tors, quark bilinears and four-quark operators. The gluonic operators can be omitted from 
consideration. This is because they transform as singlets under chiral transformations in 
the valence and sea sector, and their contribution in the effective theory will be identical 
to those from two-quark and four-quark operators of dimension-6 that do not break the 
S77(15|3)x ® SU(15\3)r chiral symmetry. We decompose the dimension-6 operators into 
three classes as follows 

W6) _ AG) , W6) , AG) (7] 

The first term consists of operators formed from only valence quark fields. When the 
valence quark masses are zero, these operators have an SU (3 1 3) SU (3 1 3) r chiral symmetry 
which is reduced to a vector symmetry when the valence masses are turned on. There is 
also a term in which breaks the SO (4) rotational invariance of Euclidean space down 

to the hypercubic group SW4. The dimension-6 terms in £^ al for unquenched theories have 
been detailed long ago lit! , for a recent discussion that addresses Ginsparg- Wilson valence 
quarks specifically, see |47[. 

The second term del consists of operators formed from only sea quark fields. The case of 



one staggered flavor (corresponding to four tastes) was considered in [18, 48]. More recently 



3 We use the super-symmetric formulation of partially quenched theories that stems back to |4flj . One could 
equivalently use the replica method 0] • 
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these results were generalized to the case of multiple flavors in |19l |21| . Terms that do not 



break taste-symmetry have the same form as operators in These operators have an 

SU(12)l® SU (12) r chiral symmetry in the massless sea quark limit, and a vector symmetry 
away from zero mass. There is a single term which reduces the rotational symmetry of taste- 
symmetric terms down to SW4. When one considers taste-symmetry breaking terms, one 
has four-quark operators of the form 

Q{T A ® T* b )V s q] \Q(r A ® T* b )V s q] (8) 

where T A denotes a Dirac matrix that acts on the spin indices of the quark fields, and T* B 
denotes a generator of SU (4) taste that acts on the taste indices of the quark fields. We treat 
the sum over staggered flavors as implicit. The presence of explicit T* B matrices breaks the 
taste-symmetry. When the collective labels A and B are unrelated, the operator respects 
SO(4) rotational invariance. On the other hand, when the Lorentz indices contained in 
label A are contracted with those in B, both taste and rotational symmetries are broken 
and consequently the distinction between the two is blurred. 

Finally operators in are operators that involve both valence and sea quark fields. 
Because of the mixed action symmetry, all of these terms are four-quark operators consisting 
of a product of a bilinear of valence quarks with a bilinear of sea quarks. These operators 
were constructed in j35[, where it was shown that only taste-singlet sea-quark bilinears are 
present. The possible Dirac structures of these bilinears are constrained by the exact chiral 
symmetry required in the valence sector, and the axial symmetry required in the sea sector. 
Only vector and axial-vector bilinears are consistent with both symmetries. Ignoring color 
structure which has no bearing in the construction of the chiral Lagrangian, we have two 
terms in 



mix 



41 = CiL (Qi^VvQ) {Ql.VsQ) + C mix {Qi^VvQ) {Q^ 5 V S Q) ■ (9) 

These terms have a SU(3\3)i <g> S77(3|3)# (g> SU{12) L ® SU{12) R chiral symmetry. 

Finally we consider the terms in C^ 7 \ We do not consider operators that have a quark 
mass insertion because these will necessarily be higher order in our power counting. There 
are no operators that consist only of Ginsparg- Wilson valence quarks because we cannot 
build parity even, chirally symmetric operators of dimension-7. As discussed in j^H, there 
are no dimension-7 operators built solely from staggered quark fields due to the staggered 
axial symmetry. Finally any operators in that consist of both staggered quarks and 
Ginsparg- Wilson quarks must occur as a product of bilinears of each fermion type. Again 
this is because there is no symmetry relating the two sectors of the mixed-action theory. 
There are no such bilinears, however, because one cannot write down a dimension-7 oper- 
ator consisting of a chirally symmetric Ginsparg- Wilson bilinear and an axially symmetric 
staggered bilinear. Thus = 0. 



B. Mesons 

In this section, we review the construction of chiral perturbation theory in the meson 
sector of this partially-quenched mixed-action theory. Partially quenched chiral perturbation 
theory was developed through a series of papers [!]J 51, 52, 53[. The construction of 



the chiral effective theory for Ginsparg- Wilson valence quarks and staggered sea quarks 
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was carried out in 35j. We will focus only on the ingredients of this theory necessary for 
determining baryon properties at next-to-leading order. 

As commented above, in the massless and continuum limit the Symanzik action has a 
graded chiral symmetry SU(15\3)l <S> SU(15\3)r. We expect this symmetry to be sponta- 
neously broken down to 5C/(15|3)y in analogy with QCD. Thus we can build an effective 
theory of mixed action PQQCD written in terms of the pseudo-Goldstone modes that emerge 
from spontaneous chiral symmetry breaking. These modes acquire masses from the explicit 
chiral symmetry breaking introduced by the quark mass term in the continuum action. 
Additionally we shall see below that the a-dependent terms in the Symanzik action also 
generate masses for some of these modes. We shall generically call the pseudo-Goldstone 
modes mesons, and collect them in an [/(15|3) matrix exponential £ 

£ = e 2 ^H = e, (10) 

that is written in terms of the meson matrix We work almost exclusively in the quark 
basis due to the fact that mesons composed of two sea quarks will enter our next-to-leading 
order calculations only as taste singlets. 

At 0(e 2 ) in the combined lattice spacing and chiral expansion the effective Lagrangian of 
partially quenched chiral perturbation theory (PQ^PT) that describes the mesons has the 
form [35[ 

f 2 1 
C = i-str (d^d^) - A str (m q tf + mj£) + -/zg (str $) 2 + a 2 V. (11) 
8 6 

With these conventions, the pion-decay constant / = 132 MeV. The potential V contains the 
effects of dimension-6 operators in the Symanzik action. It can be decomposed into three 
terms [35| 

V = U S +U' S +U V . (12) 

Above Us and U' s are taste-symmetry breaking potentials involving single trace and two trace 
operators, respectively [18l Il9j. These potentials involve only the mesons composed of two 
sea quarks and arise from C^X in the Symanzik action. The potential Uy contains all terms 
that stem from Symanzik operators involving the valence quark fields. These operators are 
in C^} and . As shown in 35||, there is only one term in the valence potential 



U v = -a 2 C mix str (t^t^) , (13) 

where T3 = Vs — Vy. This term only acts on mesons composed of one valence quark and one 
sea quark. At this order there are no contributions from operators that act on valence- valence 
mesons as their effects lead to a 2 renormalization of the lowest order parameters |47| . 

Working to tree-level, we can determine the meson masses needed below for the calcula- 
tion of baryon observables. When one expands the Lagrangian in Eq. (JTTj) to leading order, 
one finds that mesons formed of two Ginsparg- Wilson quarks QQ' have mass 

m QQ> = Jl( m Q + m Q')i ( 14 ) 

which accordingly vanishes in the chiral limit. Mesons consisting of a staggered quark Qi, 
i.e. of flavor Q and quark taste i, and a Ginsparg- Wilson quark Q' have masses 



2 4A 16 a 2 C mix _ 

m QiQ' = -p( m Q + m Q') + J 2 • ( 15 ) 
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Notice that these masses do not depend on the quark taste. Masses of mixed mesons do not 
vanish in the chiral limit because the full chiral symmetry of C is explicitly broken by 
down to the mixed action chiral symmetry. The final mesons relevant to baryon observables 
at next-to-leading order are those with two staggered quarks in a flavor-neutral, taste-singlet 
combination. In the quark basis, their masses are given by 

9 8A 64a 2 . . . . 

m QQ = j2 m Q + —( C3 + C ^ (16) 

where C 3 and C 4 are parameters entering the potential Us- 

Notice that in Eq. (JTTjl the matrix $ is not supertraceless, and we have included the si ngle t 
mass parameter /Jo- This is a device to derive the form of the flavor-neutral propagators [53]. 
Partially quenched theories have a U(1)a anomaly which renders the singlet heavy. Thus 
the form of flavor-neutral propagators can be derived in the limit that //q — > oo, and the 
singlet [here of 5C/(15|3)y] is integrated out of the low-energy theory. For the staggered 
flavor- neutral propagators, the analysis has been carried out in including the inclusion 
of 1/4-factors corresponding to the fourth-root trick. There it was shown that different 
flavor- neutral propagators exist for the taste singlet, vector and axial- vector channels. To 
the order we work, the only flavor neutrals required will be composed of two Ginsparg- 
Wilson valence quarks. Thus as with the meson mass and decay constant in mixed action 
PQxPT HI, only the taste -singlet flavor-neutral propagator is required. For a,b = u, d, s, 



this propagator is given by 

G 5ab 1 {k 2 + ml) (k 2 + m 2 „) 

where the masses of the valence- valence mesons, m 2 a and m\ hl are given in Eq. (|14|). while 
those of the sea-sea mesons, m 2 ^ and m 2 r , are given in Eq. (fl"6|) . The mass m x is defined as 
m 2 x = | + 2m 2 r ) . The flavor-neutral propagator can be conveniently rewritten as 

@VaVb = fiabPa + 7~Lab (Pa, A, Px) , (18) 



where 



Pa — — — , Ph — — — , Px 



k 2 + m 2 a ' k 2 + m 2 b ' k 2 + m 2 x ' 



H ab (A,B,C) 



ni 



j.i 



m; 
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mz 



ni 



bb 
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ml 
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X) 



(m\ - m 2 aa ) {m\ - m 



C 



(19) 



In the limit that b — > a, we require a separate form of the flavor-neutral propagator to 
handle the double pole. The functional form, however, can be related to a derivative of the 
single pole form, namely 



H aa {A,A,C) = - 



1 



d (m 2 j - m 2 aa )(m 2 rr - m 2 J ^ ( 



m 2 x )(m 2 r - m 



Xi 



dml 



(m 2 aa - m\) 



C 



[m 



x 



mi 



(20) 



keeping in mind that A = A(m 2 aa ). 
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C. Baryons 



Having reviewed the mixed action Symanzik Lagrangian and the relevant pieces of meson 
PQxPT at finite a, we now extend mixed action PQxPT to the baryon sector. In the 
continuum limit, the flavor symmetry group of the mixed action theory is SU(15\3)y. There 
are no taste-symmetry breaking interactions present when the lattice spacing is zero. Thus 
we can construct spin-1/2 and spin-3/2 baryon flavor tensors in the usual fashion [H3,|5^, 56| . 



In general we can find the dimensionality of these multiplets in SU(M\N)y using Young 
super-tableaux 57, For the spin-1/2 multiplet B^ k , we have a 



- \(M + 1)M(M - 1) + (N + 1)N(N - 1)1 + MN(M + N) (21) 
3 - 

dimensional representation of SU(M\N)y. While the spin-3/2 multiplet T^ k furnishes a 

^M(M + 1) + ^M(M + 1)(2M + 1) + ~MiV(M + N) + -(N - 1)(N - 2) - S N)0 (22) 

dimensional representation. 

For the case at hand, SU(15\3)y PQxPT, we embed the spin-| baryons in the 
1938-dimensional super-mult iplet B^ k . The spin-| baryons are embedded in the 1086- 
dimensional super-multiplet 77p . For baryon PQ%PT at next-to-leading order (and in- 
deed at next-to- next-to-leading order as well), we shall need the states in these multiplets 
consisting of at most one sea quark or at most one ghost quark. To this end, we de- 
compose the irreducible representations of SU(15\3)y into irreducible representations of 
SU(3) va \ © SU(12) sea © S , f/(3)ghost [59]. To describe these super-algebra multiplets we refer 
to their floors, where the floor number coincides with the number of bosonic ghost quarks 
contained in states. The ground floor is synonymous with zero ghost quarks. Additionally 
we refer to the levels of the multiplet |55[ to distinguish between baryon states with differing 
numbers of sea quarks. Level A baryons do not have a sea quark, level B baryons have one 
sea quark, and so on. 

The ground floor, level A of the multiplet B 1 -* consists of baryons that transform as a 
(8,1,1) under S77(3) va i © Si7(12) sca © 5 , [/(3) g host- These are the octet baryons, and are 
embedded in the tensor B l ^ k in the standard way when all of the indices are restricted to 
1-3 [60]. The first floor, level A of B^ k transforms as a (6, 1, 3) © (3, 1, 3). These states 
have been constructed explicitly in [Hf. The ground floor, level B of the B^ k multiplet 
transforms as a (6, 12, 1) © (3, 12, 1). The states constructed in [HH can be used with 
minimal modifications. One merely must re-index so that the range of g sea is extended from 
1-3 to 1-12, or merely attach a taste index to each flavor tensor. 

The situation is the same with respect to the spin-3/2 baryon tensor T^ k . Either the 
states have already been constructed in 55l . or it is trivial to extend to the present case 



those which have not. The ground floor, level A of the spin-3/2 multiplet transforms as 
a (10, 1, 1) under SU(3) va \ © S'?7(12) sea © S , f/(3) g host, and consists of the decuplet baryons 



4 This is a useful way to proceed because the external states consist of three Ginsparg- Wilson quarks 
and have a one-to-one correspondence with the baryon interpolating fields used on the lattice. On the 
other hand, baryons formed from quark tastes are complicatedly related to the usual staggered quark 
interpolating operators used on the lattice. 



S 



which are embedded in the usual fashion. The first floor, level A of 7T' J transforms as a 

(6, 1,3), while the ground floor, level B transforms as a (6, 12, 1). Thus the total number 

of baryon states relevant for calculating loop diagrams up to next-to-next-to-leading order 

is 243, which is considerably smaller than the total number of baryons in the theory 3024. 

To 0(e 2 ), the free Lagrangian for the B l ^ k and 7? fields retains the same form as in 

I — i i i " 

quenched and partially quenched theories [55|, |60j] with the addition of new lattice-spacing 
dependent terms, 

£ = i(Bv VB) - 2a M (BBM+) - 2(3 M (BM + B) - -a M (BB) str (M+) + a 2 V B 

+i (7> ■ VT^) + A (T M T M ) + 2 7M (T^M+T^) - -a M (T M T M ) str (M+) + a 2 V T . 

(23) 

The baryon potentials and W arise from the operators in of the Symanzik La- 
grangian. For our next-to-leading order calculations, we will not require the explicit form of 
either term. The effective contribution can be deduced from symmetry considerations alone. 
In the baryon Lagrangian, the mass operator is defined by 

M+^^ + M' ( 24 ) 

Above, the parameter A ~ eA x is the mass splitting between the 1938 and 1086 in the 
chiral limit. The parenthesis notation used in Eq. f!23l) is that of and is defined so that 
the contractions of flavor indices maintain proper transformations under chiral rotations. 

The Lagrangian describing the interactions of the B^ k and T^ k with the pseudo- 
Goldstone mesons is 

C = 2a (BS,BA,)+2(3 (BS^B) -2H [T V S^T V ) + J^C [{T U A U B) + {BA V T V )] . (25) 

The axial- vector and vector meson fields A^ and are defined by: A^ — | (£<9 M £t — ^d^) 
and = \ (ffyft + gg^y The latter appears in Eq. ()23)1 for the covariant derivatives of 
Bijk and %j k that both have the form 

(iy% fe = d^B ijk + (VjuB ljk + (-r^ + ^\V^) jm B imk + (-)(^)(*+*)(v (t ) tol B yB . (26) 



The vector S M is the covariant spin operator [6JJ, |62| . The interaction Lagrangian in Eq. (|23|) 
also receives lattice spacing corrections. In calculating baryon observables, however, these 
lead to effects that are 0(e 2 ) higher than the next-to-leading order results. 

The parameters that appear in the mixed-action PQ%PT Lagrangian can be related to 
those in ;\PT by matching. One realizes that QCD is contained in the fourth-root of the 
sea-sector of the theory. Thus we can relate the parameters of x?T by matching onto terms 
in the above Lagrangian restricted to one taste for each flavor of staggered sea quark. For 
definiteness, we restrict the indices to 4,8,12 corresponding to jx, h, and r x . This allows 
the usual identifications: a = |D + 2F, (3 = — ~D + F, and the remaining parameters: «m, 
Pm, &m, 1m, &m, C, and TC, all have the same numerical values as in x?T. 
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FIG. 1: Loop diagrams contributing to the octet baryon magnetic moments at 0(e). The photon 
is pictured as a wiggly line, mesons are denoted by a dashed line, and a thin solid line denotes a 
1938 baryon, while the double line denotes a 1086 baryon. 



III. BARYON MAGNETIC MOMENTS 



In this Section, we calculate the octet baryon magnetic moments in mixed action PQ%PT. 
We choose to present the calculation of magnetic moments first due to their simplicit y. R ecall 
that in the continuum limit, the baryon magnetic moments have the behavior 63[ 6^, 6o| 



/i~/io + a^+... (27) 



in the chiral expansion. Terms denoted by . . . scale with a higher power of m q , and are 
first encountered from one-loop graphs at next-to-next-to-leading order. At next-to-leading 
order then, the magnetic moments receive 0(e) contributions from loops. 

Moving away from the continuum, we are forced to address the corrections from 0(a 2 ) 
operators. While there are new a 2 -dependent terms in the baryon Lagrangian, there are 
still further terms because the vector-current operator also receives 0(a 2 ) corrections in the 
effective theory. Potential contributions from any such terms, however, scale as 0(e 2 ) and 
are relevant only at next-to-next-to-leading order. We can safely ignore such contributions. 
This situation is similar to the calculation of magnetic moments with mixed actions of 
Wilson and Ginsparg- Wilson quarks [f^, see also [i^- Thus the lattice-spacing artifacts at 
next-to-leading order can only enter through loop effects. The one-loop graphs that give a 
contribution to the octet magnetic moments are depicted in Figure ^ 

To calculate these diagrams and the tree-level contributions, we must extend the electric 
charge matrix Q, which is not uniquely denned in partially quenched theories (Hif. By 
imposing the charge matrix Q to be supertraceless in SU(15\3), no new operators involving 
the singlet component are introduced. There are various ways to accomplish this, however, 
many are not practical to implement time-wise in current lattice computations. We require 
that ghost charges equal their valence counterparts so that operator self-contractions from 
the valence quarks are completely canceled by ghosts Thus our form of the SU(15\3) 
charge matrix is 

Q = diag(g u , q d , q 8 , qj £/, q t f 7 , q r f 7 , q u , q d , q„), (28) 

where to maintain supertracelessness qj + qi + q r = 0. Taking the fourth-root of the determi- 
nant in the continuum, and making the sea quarks degenerate with the valence quarks, QCD 
is recovered only for the specific choice of charges: q u — qj — g, and qd = q s = qi = q r = — §• 
Using unphysical charges for the valence and sea quarks provides a means to access low- 
energy constants, and thereby determine the physical magnetic moments, for example. No- 
tice the choice, qj = qi = q r = 0, is allowed by supertracelessness, and would considerably 
free up computation time by eliminating all operator self-contractions. 

In PQxPT the leading contribution to the octet magnetic moments arises from two 
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dimension-5 operators |5£ 

C=—-[n a (B[5 M , S V \BQ) + tip {BIS,, S V }QB)] F, Vl (29) 



which can be matched onto the Lagrangian upon restricting the baryon field indices 
to 4, 8, 12, and taking the charges to be their physical values. In terms of the SU(3) matrix 
B containing the octet baryons, the corresponding %PT Lagrangian is 

ip 

C = ^r M r 3 "]{Q> 5 » + ^ tr ( B \. 3 *> s "1 t Q ' S D (3°) 

and hence we find \i a = \\i>d + 2/^?, and /^g = — + A*f by matching. These operators 
contribute to the magnetic moments at O(e ). There are additional operators that make 
contributions of 0(e) that have identical flavor structure, and differ only by the insertion 
of A/A x . These operators come with new low-energy constants and are allowed because 
the mass-splitting parameter is a chiral singlet. These A-dependent operators also function 
as counter-terms for the one-loop divergences. We shall not keep these operators explicitly, 
and merely treat the leading low-energy constants as arbitrary linear functions of A. The 
ability to determine coefficients for the A-dependent operators requires the ability to vary 
A, and for this reason we treat the A dependence only implicitly. 

Having spelled out the operators for tree-level contributions to the magnetic moments, 
we now comment on the computation of the loop diagrams. Because flavor-neutral mesons 
remain charge neutral, even with arbitrary quark charges, there are no hairpin contributions 
at next-to-leading order, see Figure [TJ Loop graphs containing intermediate state baryons 
with all valence quarks, or two valence quarks and one ghost quark are identical to those 
calculated for SU(6\3) in jH^] (modulo the different choice of electric charge matrix) because 
the ground floor, level A, and first floor, level A transform identically in each theory (both 
are singlets in the sea sector where the group structure is different). Finally there are the 
diagrams with baryons consisting of two valence quarks and one sea quark. The correspond- 
ing loop mesons are mixed mesons; they have one valence quark and one sea quark. The 
propagation of these modes does not break taste-symmetry. This is because the Lagrangian 
in the mixed sector does not contain explicit taste matrices only the matrix T3. Thus the 
mixed meson propagator obeys |35| 

*^,^) = (3D 

where Q iy Qj label the staggered quark flavors and tastes, and Q', Q' label the Ginsparg- 
Wilson quark flavors. The mixed meson mass mn-Qi is given in Eq. (|15jl . and is independent 
of the quark taste i. Now because there is no taste changing in the loop, and the relevant 
meson and baryon masses are quark taste independent, we obtain a factor of four for each 
staggered quark flavor that can propagate in the loop. This factor is then canceled by the 
fourth-root trick. Explicit calculation verifies this. The baryon magnetic moments hence 
have the same form as in SU (6\3) PQxPT (Eif, and depend only on one new low-energy 
constant C m i x . This low-energy constant appears in the loop meson masses that contain a 
sea quark and a valence quark. 



5 Here we use — d,A u — d v A p 
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TABLE I: Tree-level coefficients for the magnetic moments of octet baryons. 



P 
n 

£+ 
S° 

s- 

H° 
A 

AS 



Q 



q u 



q s 
q a 



q u - qd 
o 

-q u + qd 
-qd + q s 
-q u + q s 

o 

o 



an 



q u + q s 

qd + q s 

q u + qd 
q u + qd 
q u + qd 
qd + q s 
q u + q s 
\{q u + qd + ^q s 



7s 



Combining the tree-level and one-loop graphs, we have 



T{m^ A, \i) + - 



(32) 



where we have carried out the spin algebra in D-dimensions. 
appearing in the above expression is defined by 



The non-analytic function 



jF(m, 5, /i) = VS 2 — m 2 log 



5 — V<5 2 + m 2 + it 



m 



5 + V6 2 + m 2 + ie 



<51og— . 
fi 2 



(33) 



The coefficients Q, and otp for the tree- level diagrams are listed in Table HJ Ordinarily Q 
is the baryon charge; this is no longer the case with the quark electric charge matrix Q 
given in Eq. (|2Sjl. The computed values for the (3^ and /% coefficients are listed for the 
octet baryons in Tables ITT MlXl The corresponding values of these coefficients for the AS 
transition moment are given in Table IXl In each table we have listed the values corresponding 
to loop mesons that has mass m^. In these diagrams there are valence- valence mesons, with 
masses given in Eq. (|14jh and valence-sea mesons, with masses given in Eq. (|15p. If a 
particular meson is not listed then the values for /3^, and (3L are zero. 

The charges qj and qi do not appear explicitly in these tables because, in the isospin limit, 
they always come in the combination qj + qi and this is identical to — q r . Again we remark 
that the supertracelessness of Q is maintained by the computation-time simplifying choice 
qj = qi = q r = 0. For this choice of charges, one is not computing the physical magnetic 
moments; however, the unphysical moments determined are sensitive to the physical low- 
energy constants, hence predictions can be made. Additionally one has complete freedom 
to adjust the charges q u , qd, and q s as these do not contribute to the supertrace of Q. One 
can, for example, choose these charges to isolate the low-energy constants at tree-level, see 
Table QJ or to simplify the chiral extrapolation by eliminating particular loop mesons. 

The form of the continuum extrapolation of baryon magnetic moments is thus highly 
constrained, especially since the only new low-energy constant at this order, C m i x , could be 
determined independently from the masses of mixed mesons [35] . 
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TABLE II: The coefficients fix, and (3' x for the proton. 




h ft 


7T 

ju 
ru 


-^D 2 {q u -q d ) -\{q u -qd) 
-l(D 2 + 3F 2 )q u - 2(D - F) 2 q d - (f D 2 - 2DF + 3F 2 ) q r ±{2q u + 4q d + 3q r ) 
-l(D 2 + 3F 2 )q u - (D - F) 2 q d + (\D 2 - 2DF + 3F 2 ) q r ±(q u + 2q d - 3q r ) 


TABLE III: The coefficients fy, and ^ for the neutron. 






7T 

ru 


±D 2 (q u -q d ) \{qu-qd) 
-2(F) - F) 2 q u ~I(D 2 + 3F 2 )q d - (§D 2 - 2DF + 3F 2 ) q r ±(4q u + 2q d + 3q r ) 
-(D - F) 2 q u -l(D 2 + 3F 2 )q d + (§D 2 - 2DF + 3F 2 ) q r ±(2q u + q d - 3q r ) 


TABLE IV: The coefficients and ^ for the £+. 




P<t> 


if 

TO 

is 

rs 


-^D 2 {q u -q s ) -\{qu-q s ) 
-§(D 2 + 3F 2 ){2q u +q r ) ^(2q u + q r ) 
-l(D 2 + 3F 2 ){q u -q r ) ±{q u - q r ) 
-{D - F) 2 {2q s +q r ) \{2q s + q r ) 
-(D-F) 2 (q s -q r ) ±(g s - q r ) 


TABLE V: The coefficients /fy, and for the E°. 




h P* 


if 

ju 

TO 

is 

rs 


-\D 2 (q u + q d - 2q s ) - -^(q u + q d - 2q s ) 
~1{D 2 + 3F 2 ){q u + q d + q r ) ±(q u + q d + q r ) 
-\{D 2 + 3F 2 )(q u + q d - 2q r ) ±(q u + q d - 2q r ) 
-(D-F) 2 (2q s + q r ) ±(2q s + q r ) 
-(D-F) 2 (q s -q r ) \{q s - q r ) 


TABLE VI: The coefficients /fy, and for the E~. 






if 

ju 
ru 

is 

rs 


-^D 2 (q d -q s ) -\{qd-Qs) 
-l(D 2 + 3F 2 )(2q d + q r ) Tg{2q d + q r ) 
-l(D 2 + 3F 2 )(q d -q r ) ig(?d - <&■) 
-{D-F) 2 (2q s + q r ) l{2q s + q r ) 
-(D-F) 2 (q s -q r ) ±(q s - q r ) 
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TABLE VII: The coefficients 0+, and for the S°. 






P<p 




K 

ju 
ru 
js 
rs 


^D 2 (q u - q s ) 
-(D-F) 2 (2q u +q r ) 
-(£> - F) 2 (q u - q r ) 
-§(£> 2 + 3F 2 )(2q s +q r ) 
-l(D 2 +3F 2 )(q s -q r ) 


\(q u - q s ) 

\(2q u + q r ) 

Ulu - Qr) 

j g {2q s + q r ) 



TABLE VIII: The coefficients and ft, for the H". 



</> 


h ft. 


K 

ju 
ru 
js 
rs 


lD 2 (q d -q s ) 
-(D-F) 2 (2q d + q r ) 
-(D-F) 2 (q d -q r ) 
-l(D 2 + 3F 2 )(2q s +q r ) 
-l(D 2 +3F 2 )(q s -q r ) 


Hid ~ Qs) 
l(2q d + q r ) 

\{<ld - q r ) 
jg{2q s + q r ) 

T§(qs-q r ) 



IV. BARYON MASSES 



In this Section we determine the octet baryon masses in mixed action PQ^PT. Let us 
first recall the behavior of the baryon mass in the continuum, and near the chiral limit [l^, 

Mb ~ M + a m q + (3 m 3 J 2 + . . . , (34) 

where we have retained the leading non-analytic piece and the . . . denotes terms with higher 
powers of the quark mass. The non-analytic term above stems from the one-loop diagrams 
shown in Figure El A mass calculation at next-to-leading order is 0(e 3 ) in our power 
counting. To perform a complete 0(e 3 ) calculation, we must retain the a 2 -dependent piece 
of the loop meson masses and evaluate the potential Vb in Eq. f!23D at tree level. Remember 
that there are no a 3 terms in the Lagrangian. 

Without writing down the explicit form of the baryon potential, we can deduce the net 
contribution from Vb to the octet masses at tree level. Because the external states all 
involve Ginsparg- Wilson quarks, the action of the potential is necessarily a taste-singlet at 



TABLE IX: The coefficients fy, and ft, for the A. 









K 

ju 
ru 
js 
rs 


lD 2 {q u + q d -2q s ) 
-l(7D 2 - 12DF + 9F 2 )(q u + q d + q r ) 
-l(7D 2 - YIDF + 9F 2 )(q u + q d - 2q r ) 
-l{D 2 + 6DF + 9F 2 ){2q s + q r ) 
-\(D 2 +QDF + 9F 2 )(q s -q r ) 


j^iqu + q d - 2q s ) 
\{q u + q d + q r ) 
j2(q u + q d - 2<2v) 
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TABLE X: The coefficients 0^, and for the AS transition. 








7T 

K 

ju 
ru 


-s73^ 2 (9«-9d) 

^j(D 2 — 2>DF)(q u — q d ) 
^(D 2 - 3DF)(q u - q d ) 


12V z(qu q d ) 




FIG. 2: Loop diagrams contributing to the octet baryon masses at C(e 3 ). Mesons are denoted by 
a dashed line, flavor neutrals (hairpins) by a crossed dashed line, and a thin solid line denotes a 
1938 baryon, while the double line denotes a 1086 baryon. 

tree level. Furthermore since insertion of the mass matrix is at 0(e A ), there are no valence 
flavor matrices around. This implies that the evaluation of the baryon potential is identical 
for all octet states: in essence all Ginsparg- Wilson quarks are identical in the potential at this 
order. While there are terms in the potential that violate 5*0(4) rotational invariance |74| . 
the mass only picks up an indiscernible additive shift from such terms [it is the dispersion 
relation that is sensitive to 50(4) breaking]. Therefore the effect of the potential at 0(a 2 ) 
can be described by just one operator 

V B = -C (BB) . (35) 

This term and the familiar tree-level terms that are linear in the quark masses are straightfor- 
wardly evaluated. We must keep in mind that qm, Pm, and % are to be treated as arbitrary 
linear functions of A/A x . The linear terms, for example, function in part as counter-terms 
for the one-loop diagrams. 

Now we turn to the loop diagrams in Figure El As with the computation of the magnetic 
moments, the diagrams involving baryons with three valence quar ks, or two valence quarks 
and one ghost quark are identical to those in 577(6 1 3) PQxPT |55j, because these particular 
floors and levels of each multiplet are identical. Furthermore, diagrams with one sea quark 
flowing in the loop are evaluated in precisely the same way as before. The mixed-meson 
propagators are diagonal in taste and flavor, while the meson and baryon masses are taste- 
independent. Consequently the loops come with a four-fold degeneracy when one sums over 
the sea-quark tastes. This degeneracy is canceled by the factor of 1/4 to implement the 
fourth-root trick. The only new contribution as compared to the magnetic moments is that 
of flavor-neutral modes. But at this order, the flavor-neutral mesons which propagate in 
loops are those consisting of two Ginsparg- Wilson valence quarks. These are necessarily 
taste-singlets; their propagators appear in Eq. (fTTj). and involve taste-singlet sea-sea meson 
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TABLE XI: Tree-level coefficients for the octet baryon masses. 





C u 


c d 


N 


OCM + 0M 





£ 


\{ba M + 2f3 M ) 


\(a M + 4/?m) 


A 


\{a M + WM) 


\olm 




\(olm + 4(3 m) 


l(5a M + 2(3 M ) 



masses, which are a 2 -dependent. 

Assembling the tree-level results with the loop-diagrams, we find 

Mb = M — 2m u C u — 2m s C s — 2crM(2mj + m r ) — a 2 Co 



1 



W 2 
C 2 



A <t> m \ + A W'M?{m^ 



8tt 2 / 2 



22 B <j) F{m^ A, /x) +2j%F(m^,m^,A,//) 



(36) 



The non-analytic functions appearing in the expression for the octet baryon masses are 
defined by 



F(m,5,jj) = [m 2 —5 2 ) V5 2 — m 2 lo{ 



5 - V5 2 + m 2 + 



IE 



5 + y/ 8 2 + m 2 + it 



5 log — )--Sm 2 log—, (37) 



and 



F{m^ mp, 5, /x) = H^F (m^ 5, /x), F(m^>, 6, /x), F(m x , 6, /x)], 



(38) 
(39) 



for the flavor-neutral contributions. 

In Table 1X11 we list the coefficients C u and C s for the octet baryons. In Tables 1X111 and 
1X1111 we list the loop coefficients A^, A^/, and B^i. These coefficients are grouped 
according to loop mesons with mass m^, and for flavor-neutral contributions are grouped 
according to pairs of quark-basis flavor-neutral mesons. Notice the sum on 00' runs over 
VuVu, VuVs, an d rjsfjs to avoid double counting. 

The lattice spacing dependence of the octet masses too is highly constrained: there 
are only three free parameters: Co, which is the same for all octet baryons; C m ; x , which 
alternately can be determined from mixed meson masses; and the combination of parameters 
C 3 + C 4 , which is already constrained from staggered meson lattice data p. 



V. SUMMARY 



Above we have included baryons into mixed-action partially-quenched chiral perturbation 
theory for Ginsparg- Wilson valence quarks and staggered sea quarks. Working at next-to- 
leading order, we determined the lattice spacing artifacts for the octet baryon magnetic 
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TABLE XII: The coefficients and A^i in PQxPT. Coefficients are listed for the baryon octet, 
and for A§ are grouped into contributions from loop mesons with mass m^, while for A^> are 
grouped into contributions from pairs of quark-basis rj q mesons. 







A^ 








7T 


K 


Vs 




N 


-l(D 2 -3DF) 










E 


-*(D 2 -3F 2 ) 


-|(,D 2 - 6DF + 3F 2 ) 







A 


— 2 (D 2 — ^'2,DF + 9F 2 \ 


_2(cr)2 _cnr_ ni?2\ 


n 









-§(L> 2 - 6DF + 3F 2 ) 


_| (jD 2_ 3f2) 










js 


rs 


N 


§(5£> 2 - 6DF + 9F 2 ) 


±(5L> 2 - 6DF + 9F 2 ) 








E 


l(D 2 + 3F 2 ) 


l(D 2 + 3F 2 ) 


2{D-F) 2 


(D - F) 2 


A 


|(7 D 2 - 12DF + 9F 2 ) 


§(7£> 2 - 12DF + 9F 2 ) 


l(D + 3F) 2 


\{D + 3F) 2 




2{D - F) 2 


{D - F) 2 


l(D 2 + 3F 2 ) 


l(D 2 + 3F 2 ) 






A^i 








VuVu 




VsVs 




N 


(D - 3F) 2 










E 


4F 2 ' 


-4(DF- F 2 ) 


{D - F) 2 




A 


l(2D-3F) 2 


-§(2£> 2 +3DF-9F 2 ) 


\(D + 3F) 2 






(D - F) 2 


-4(DF - F 2 ) 


4i?2 





TABLE XIII: The coefficients and B^i in PQxPT. Coefficients are listed for the octet baryons, 
and for are grouped into contributions from loop mesons with mass m^, while for B^> are 
grouped into contributions from pairs of quark-basis rj q mesons. 
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Vs 


ju 


ru 


js 


rs 
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3 








2 
3 
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2 
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9 


9 
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moments and masses. The recipe for adding finite volume corrections is discussed in the 
Appendix. 

To 0(e 3 ), baryon masses depend on the lattice spacing via three parameters. The first 
is Co, which is a representative coefficient of local a 2 operators in the baryon Lagrangian. 
The contribution from C is the same for all members of the baryon octet. The second 
parameter is C m i x , which affects the masses of mesons made from a Ginsparg- Wilson quark 
and a staggered quark. The third parameter is a combination of low-energy constants, 
C3+C4, which governs the mass of taste-singlet staggered mesons. To 0(e), baryon magnetic 
moments only depend on the lattice spacing through loop-meson masses. Furthermore there 
is only one free parameter involved at next-to- leading order, C mix . 
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A simplifying feature of the mixed action baryon theory is that at next-to-leading order, 
there is no taste-symmetry violation. Beyond next-to-leading order, however, this is no 
longer the case. Taste-symmetry breaking interactions start at 0(e A ) in our power count- 
ing. They arise from two sources. Contributions that scale as 0{m 2 q \ogm q ) can lead to 
taste-symmetry violation. Expanding out the sigma terms in the baryon Lagrangian (|23|) 
leads to contributions to the baryon masses that scale as m 2 logm g , and involve completely 
disconnected meson loops. In PQ^PT these loops involve only mesons formed from two 
sea-quarks |75|, [Til [77j . For staggered sea quarks this becomes a sum over the various me- 
son tastes (vector, axial- vector,...) and no longer is there a cancellation of factors of 1/4 
inserted from the fourth-root trick. Explicit taste-symmetry violation occurs from operators 
that scale as 0(a 2 m q log m q ). These arise from the loop-diagrams generated by terms in 
Vb (for which one would need to decompose Cq into various contributing terms). Terms in 
Vb containing taste-spurions will generate taste-symmetry violation at the one-loop level. 
These are the only possible sources for taste-symmetry violation at next-to-next-to-leading 
order. The remainder of contributions to the baryon masses and other observables arise 
from valence-valence or valence-sea mesons in the absence of explicit taste matrices. These 
do not violate taste-symmetry as we demonstrated above. 

Nonetheless the issue of addressing lattice spacing corrections to baryon observables cal- 
culated in mixed-action lattice QCD is very tractable. The chiral symmetry properties of 
Ginsparg- Wilson valence quarks effectively suppress the taste-symmetry violation from the 
staggered sea. While taste-symmetry violation does occur at next-to-next-to-leading order, 
expressions at leading-order are taste-symmetric and involve only a few new parameters. 
Mixed- act ion simulations in the baryon sector are a timely way of getting physical observ- 
ables from lattice QCD. 
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FINITE VOLUME CORRECTIONS 



In this Appendix, we assemble the relevant formulae for finite volume corrections. This 
is a trivial extension of Ref. and is included for completeness. Let L denote the size of 



the cubic box in one spatial direction. We assume the lattice simulations are carried out in 
a region of parameter space where chiral physics lives inside the box, i.e. for fL ^> 1, and 
further that we are in the p- regime of chiral perturbation theory, m^L > 1. In this case, the 
Poisson formula can be used to cast mode sums from loop diagrams into the infinite volume 
results plus finite volume modifications. We list these modifications for the octet baryon 
masses and magnetic moments. For finite volume corrections in mixed action simulations, 
we imagine that the Ginsparg- Wilson valence quarks will be the lightest. Thus modulo 
possible cancellations from additive a 2 mass shifts for the valence-sea and sea-sea mesons, 
the valence-valence pions should dominate the finite volume corrections. 
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The finite volume corrections to the magnetic moments have the form |7j 

M B 



6tt 2 / 2 



^[^K,O)+/3;C 2 y(m ,A)] 



where 



y{m,A) = / dX^ hK (p A \n\L) - (3 A \n\L K^^L) 



n-0 



(40) 



(41) 



with (3\ = m 2 + 2 AX + X 2 and the K n (x) are modified Bessel functions. The coefficients (3$ 
and (3'^ are listed for the octet baryons in Tables ITTMlXI The coefficients for the AE°transition 
moment appear in Table IXl 

The finite volume corrections to the masses have the form 176 



5 L M B 



4tt/ 2 
C 2 

47T 2 / 2 



y^^/C(w^,0) + A^Kjm^ m# , 0) 
y^-B^/C(m^,, A) + B^Kjm^ 7ry , A) 



where 



and 



AT(m, A) 



^l(W) 

(3 A \n\L 



K (f3 A \n\L) 



K.(m^, mp, A) = [/C(m , A), /C(m^, A), K,(m x , A)] . 



(42) 



(43) 



(44) 



The coefficients and A^< are listed for the octet baryons in Table IXII1 while and 
Bstt,' appear in Table IXIII1 
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